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In  order  to  aeseaa  the  potential  of  multi  stage  pm- launched 
rockets,  a  stud?  of  tkS  vehicle  and  trajectory  p*r*a*t«'  was 
undertaken.  A  digital  computer  program  for  trajectories  was 

t 

written  end  was  used  in  an  experimental  wanner  to  approach  optimum 
performance  within  various  sets  of  restricting  assimptions.  The 
approach  was  found  to  be  effective  and  a  useful  orbital  potential 
urns  deswnstrated  with  reasonable  design  parameters. 
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INTRODUCTION 

During  recent  years  interest  has  been  growing  in  the  use  of 
guns  to  deliver  payloads  to  high  altitudes.  A  substantial  amount 
of  development  in  this  area  has  been  carried  out  by  the  Ballistic. 
Reoearch  Laboratories  of  the  U.S.  Army  and  by  McGill  University 
with  supDcrt.  from  the  U.S.  Army  through  the  Army  Research  Office 
and  B.R.L.  (Refs.  1  &  2).  The  latter  effort,  has  been  concerned 
with  a  16  inch  gun  located  on  Barbados  (Refs.  3  &  4).  The  gun  is 
smooth  bored  and  is  mounted  to  allow  firing  elevations  from  about 
60°  to  near  vertical.  Initial  developments  were  concerned  with 
sub-caliber  fin  stabilized  saboted  ballistic  projectiles.  This 
class  cf  vehicle  has  been  fairly  well  proven  at  this  stage,  and 
operations  are  phasing  into  scientific  data  gathering  flights. 

At  the  same  time,  a  sub-caliber  rocket  assisted  vehicle  is  being 
developed.  (Ref.  5)  Ibis  vehicle  will  eventually  be  useful  for 
research  purposes  with  higher  performance  requirements  than  can  be 
achieved  with  pure  ballistic  vehicles.  It  is  also  a  stage  in  the 
development  of  a  multi-stage  full  bore  rocket  powered  vehicle 
with  orbital  capabilities. 

In  order  to  provide  a  basis  for  orientation  and  design 
guidance  of  a  program  to  develop  such  a  vehicle,  it  was  necessary 
to  indicate  achievable  performance,  based  on  suitable  assumptions 
and  corresponding  design  requirements. 


In  particular  it  was  naceaaary  to  determine  the  type  of 
trajectory  which  ahould  be  flown  and  the  corresponding  citing 
of  the  rocket  stages.  It  waa  decided  that  experimental  use  of 
a  computer  program  for  trajectory  calculation  waa  a  promiaing 
approach.  The  intention  waa  to  sub-optimise  what  appeared  to 
be  the  more  important  parametera,  and  then  to  determine  the 
effects  of  other  variables,  with  aome  checking  to  enaure  that  no 
important  croas  effects  were  neglected. 

Although  many  computer  programa  for  trajectory  computation 
are  available,  all  have  certain  restriction*  and  many  Include 
detail  which  would  not  be  required  in  a  study  of  this  nature. 

It  was  decided  that  a  program  designed  specifically  for  this 
task  would  provide  benefits  far  out  weighing  the  time  end 
expense  of  writing  and  debugging  the  program,  particularly  alnce 
a  large  high  speed  computer  (IBM  7040)  with  an  advanced  and 
efficient  language  system  (Fortran  IV)  was  available  at  McGill 
at  reasonable  coat.  Accordingly,  the  progrsm  waa  written  and 
a  full  description  of  it  is  Included  in  this  report,  It  has 
provan  to  be  extremely  useful  for  this  study,  and  it  can  be 
modified  with  little  effort  for  other  uses,  such  as  the 
production  of  tracking  data. 

The  technique  of  experimental  optimisation  has  proven  to  be 
quite  effective.  Usually  the  combined  effects  of  two  parameters 
could  be  determined  with  sufficient  accuracy  for  these  studies 
from  two  computer  runs  comprising  nins  to  twelve  cases 
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(tra Jeotories  and  vehicle*)  each.  This  would  c;st  of  the  order 
of  $20.  In  computer,  tire  ar.d  about.  the  *ame  fee  preparation, 
and  analysis.  All  twld  about.  500  ca*e*  have  been  run  to  data 
and  all  of  the  variable*  considered  to  be  of  major  importance 
have  been  Investigated.  Although  it  is  not  claimed  that 
performance  increase*  could  rot;  be  achieved  by  further 
manipulations  of  the  parameters,  ir.  is  felt.  that,  the  results  to 
date  represent  a  close  approximation  to  the  ultimate,  such  that 
any  improvement  in  performance  could  be  considered  a  detail 
refinement;  and  not  a  fundamental  change  in  design,  philosophy. 
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DESCRIPTION  OF  COMPUTER  PROGRAM 

The  program  calculates  trajectories  for  multi-stage  gun 
launched  rockets.  It  is  restricted  to  the  three  translational 
degrees  of  freedom,  ar.d  to  axially  symmetric  vehicles  with 
constant  thrust  motors.  The  earth  ia  assumed  to  be  spherical, 
and  the  following  functionel  relationships  are  pre-programmed: 
saro  lift  drag  coefficient  end  normal  force  coefficient-angle 
of  attack  derivative  vs  Mach  number  (Fig.  1),  ambient  pressure 
and  sonic  speed  ve  altitude  (Figs.  2  &  3).  Internal  ballistics 
are  not  calculated.  Ignition,  of  each  stage  occurs  at  a 
requested  value  of  altitude,  absolute  flight  path  angle,  or 
time  after  burnout  of  the  previous  stage  (or  launch,  in  the 
case  of  the  first  stage).  The  attitude  of  each  stage  during 
engine  firing  may  be  in  the  direction  of  the  relative  velocity 
(as  for  an  aerodynamical ly  stabilized  vehicle),  in  the  direction 
of  the  absolute  velocity  (tangent  to  the  absolute  trajectory),  a 
fixed  attitude  tangent  to  the  absolute  trajectory  at  ignition,  or 
a  fixed  attitude  specified  by  elevation  and  azimuth  at  ignition. 

The  equations  of  motion  ere  written  with  respect  to  axes 
defined  as  follows: 

U  -  along  the  radius  vector 

V  -  along  the  horizontal  velocity  vector 

W  -  to  fora  a  right,  handed  orthogonal  triad 

Position  is  related  through  direction  cosines  to  a  set  of  non- 
rotating  axes  defined  as  follows: 

X  -  in  the  equatorial  plane  toward  the  launch 
longitude  at  the  time  of  launch 

Y  -  to  form  a  right  handed  orthogonal  triad 

Z  -  normal  to  the  equatorial  plane  toward  the 
north 
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The  equation*  are  Integrated  by  the  Kutta-Gill  method  (Ref.  6). 

The  symbol*  used  are  defined  in  Appendix  I.  The  mathematical 
relationship*  are  listed  in  Appendix  II,  corresponding  in 
order  to  the  complete  program  listing  in  Appendix  III. 

The  inputs  t:o  the  program  are  listed  in  Appendix  IV,  along 
with  the  card  formats,  and  the  instructions  for  deck  make-up  are 
in  Appendix  V.  There  are  twelve  general  inputs  on  the  first  data 
card  for  each  case,  followed  by  a  card  for  each  rocket  stage 
containing  eleven  inputs.  All  data  fields  are  five  columns  and 
decimal  points  need  not  be  punched. 

The  outputs  of  the  pr-gram  are  listed  in  Appendix  VI,  which 
also  includes  a  sample  print-out;.  Inputs  are  listed  automatically 
at  the  start  of  each  case.  In  addition,  messages  are  provided 
at  ignition,  burnout,  apogee  and  impact.  Latitude  and  longitude 
are  also  provided  at  ignition,  burnout,  and  impact;.  Perigee  and 
apogee  heights  are  calculated  if  burnout  occurs  at  an  altitude 
greater  than  250,000  feet  and  the  orbit  is  net  hyperbolic.  The 
program  continues  pest  final  stage  burnout  to  Impact  only  if  the 
calculated  perigee  height  at.  final  stage  burnout  is  below  50  ran. 
Each  case  is  terminated  by  a  listing  of  the  errors  in  the  U  and 
W  unit  vectors,  and  the  value  of  the  UW  unit:  vector  scalar 
product,  followed  by  messages  concerning  calculation  error 
conditions  (underflow  may  be  neglected)  and  computing  time. 
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EXPERIMENTS  AND  RESULTS 

The  work  completed  so  far  has  been  restricted  to  orbital 
c«p«|fXViUUL!  UX.  a : Uxraas  taga  rocke t ,  with' variations  in  gun 
•lavatlon,  first  and  second  stags  ignition  points,  first  and 
second  stags  thrust  direction,  and  third  stags  weight  (with 
fixed  all  up  vel£it).  It  is  obvious  that  specific  impulse, 
mass  fraction  and  mussla  velocity  are  major  factors  in  determining 
performance,  but  their  primary  contribution  is  to  the  ideal 
velocity  Increment.  In  fact,  specific  impulse  and  mass  fraction 
may  be  traded  off  against  one  another  or  against  payload  with 
cooperative  impunity  by  maintaining  the  ideal  velocity  increment 
of  each  stage  at  the  original  value.  Errors  would  arise  due  to 
relatively  slight  changes  in  ballistic  coefficient  after  first 
stage  ignition,  but  most  of  the  drag  loss  has  been  incurred 
prior  to  this  and  these  errors  should  not  be  significant. 

Similar  calculations  with  muzzle  velocity  as  an  additional 
parameter  are  not  likely  to  be  as  accurate  due  to  drag  losses  at 
low  altitude  and  trajectory  changes,  but  could  be  attempted  as  a 
first  approximation  (i.e.  an  increase  in  muzzle  velocity  could 
be  interpreted  as  a  decrease  in  all  up  weight,  an  increase  in 
payload  weight,  or  decreases  in  specific  impulse  or  mass  fraction, 
by  maintaining  the  same  total  ideal  velocity  increment). 

The  general  approach  consisted  of  selecting  arbitrarily  an 
apparently  reasonable  set  of  parameters,  and  then  varying  two 
parameters  at  a  time  with  three  or  four  values  each. 
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One  of  the  variables  was  always  gun  elevation,  and  the  apogee 
and  perigee  heights  after  third  stage  burnout  were  plotted 
against  gun  elevation  to  determine  the  achievable  circular 
orbit  heights  and  the  corresponding  gun  elevations.  Typical 
remits  are  sketched  in  Fig.  4.  Apogee  and  perigee  heists 
are  represented  by  two  curves  against  gun  elevation.  Where 
circular  orbits  are  possible,  the  curves  cross,  and  exchange 
identification  because  of  the  definitions  of  apogee  and  perigee. 
Both  curves  however  are  continuous  provided  that  the  flight  path 
angle  at  third  stage  burnout  ia  sero  (i.e.  the  velocity  vector 
at  this  point  ia  horizontal).  One  curve  represents  the  hel£it 
at  burnout,  the  other  represents  the  height  half-way  around 
the  earth.  Where  the  ignition  time  and  throat  direction  of  the 
third  stage  are  not  adjusted  correctly  to  produce  zero  flight 
path  angle  et  burnout,  the  curves  cannot  quits  meet  and  are 
modified  in  the  region  near  the  intersection  auch  that  apogaa 
and  perigee  heights  respectively  are  continuous  functions  of 
gun  elevation. 

Ihe  following  parameters  were  hsld  constant  at  tha 
Indicated  values  (latitude,  longitude,  azimuth, 
launch  weight  and  cross  section  area  rafar  to  rhs  present  16.4 
inch  gun  installation  at  Barbados). 


Muzzle  Height 

Latitude 

Longitude 

Admuth 

Muzzle  Velocity 


150  ft. 

13.07  deg. 
•59.48  deg. 
118.23  deg. 
4500  ft.  sec. 


Launch  Weight 

2000  lbs. 

Nuaber  of  Stages 

3 

Cross  Section  Area 

1.478  ft. 

Nozzle  Exit  Area 

0.785  ft. 

Burning  Time 

20  sec. 

Mass  Fraction 

0.8 

Spsclflc  Impulse 

300  tec. 

Jhe  third  stage  was  constrained  to  fire  horizontally,  and 
the  absolute  flight  path  angle  at  ignition  was  adjusted  to  the 
nearest  0.01  dag.  to  keep  the  absolute  flight  path  angle  at 
burnout  within  0.01  deg. 

The  initial  set  of  parameters  included  a  fairly  ligit 
payload  of  16  lbs. ,  geometric  staging  (equal  ideal  velocity 
lncraatnts),  first  stage  i£iltion  at  100,000  ft..,  and  first  and 
second  stage  thrust  direction  type  0  (see  list  of  firing 
Indicator  values,  Appendix  IV).  The  first  parameter  varied  in 
addition  to  gun  elevation  was  the  second  stage  Ignition  point, 
and  the  results  are  shown  in  Fig.  5.  It  is  obvious  that  the 
second  stage  should  be  ignited  as  soon  as  possible  after  first 
stags  burnout.  A  similar  effect  for  first  stage  Ignition  la 
indicated  in  Fig.  6  (16  Iba.  payload)  and  Fig.  7  (40  lbs. 
payload),  although  in  the  latter  case  the  Increase  in  drag 
loaaaa  due  to  ignition  at  low  altitudes  evidently  is  a  limiting 
factor,  amphasized  by  the  low  launch  angle.  Fig.  6  also  shows 
checkpoints  indicating  the  effect  of  second  stage  ignition 
delay  and  first  stage  firing  direction.  It.  is  apparent  that 
psrformance  la  not  extremely  sensitive  “o  these  factors. 
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The  effect  of  second  stage  ignition  delay  was  rechecked 
for  a  40  lb.  payload,  and  at  the  same  time  the  effect  of 
second  stage  firing  angle  was  investigated  (Fig.  8).  Again, 
the  effects  are  not  of  major  importance  within  reasonable 
ranges  of  values;  the  desirability  of  early  second  stage 
firing  is  still  evident,  and  it  is  also  clear  that  a  g-turn 
trajectory  is  close  to  the  optimum. 

Attention  was  then  focussed  on  stage  sizing.  Since  the 
first  and  second  stages  should  be  fired  in  quick  succession, 
it  was  assumed  that  geometric  sizing  (equal  velocity 
increments)  would  be  near  optimum  for  these  stages,  and  the 
more  significant  effect  would  be  the  weight  cl  the  third 
stage  (and  thus  the  contribution  of  the  third  stage  to  the 
total  ideal  velocity  increment).  This  effect  is  shown  in 
Figs.  9,  10  &  11.  Again,  performance  is  not  extremely 
sensitive  to  parameter  variations.  The  optimized  configuration 
yields  a  slight  increase  in  circular  orbit  height,  and  exhibits 
rather  large  variations  from  geometric  staging  (Fig.  12). 
Comparison  of  Fig. 6  A  Fig. 7  indicates  that  the  first  stage 
ignition  altitude  should  increase  as  the  payload  increases. 

Thus  the  performance  indicated  in  Fig.  11  (constant  first 
stage  ignition  altitude  at  25,000  ft)  is  probably  leaa  than 
that  achievable  with  the  heavier  payloads.  In  fact,  the  point 
at  50  lb.  payload  required  lowering  the  second  stage  firing 
angle,  and  stage  weights  were  not  optimized  since  it  probably 


would  have  been  neceeaary  to  consider  the  weights  of  all 
three  stages  as  well  as  the  second  stage  firing  angle 
and  the  gun  elevation.  It  may  be  possible  to  retain  a 
g-tum  trajectory  with  later  ignition. 

It  is  interesting  to  note  that  as  the  parameter  values 

approach  the  optimum,  tr»e  gun  elevation  approaches  the 

o 

neighbourhood  of  45  ,  and  slopes  of  the  curves  of 
apogee  and  perigee  height  va.  gun  elevation  approach  equal 
magnitude  (with  opposite  signs).  Conversely,  for  any 
condition  which  la  far  from  optimum,  these  characteristics 
are  not  apparent. 

It  is  of  interest  tc  examine  the  losses  involved  in  the 
trajectories  considered.  Considering  two  payloads  (16  and 
40  lb.),  optimised  third  stage  weight  and  first  stage 
ignition  at  25,900  ft.  (not  optimum  fcr  40  lb.  payload)  we 
have  the  following  energy  balances:- 


Payload  Weight  -  Lb.  16  40 

Input 


Hassle  Velocity 

ft/sec 

4500 

4500 

locket  Velocity 

Increment  ft/sec 

29019 

25338 

Total  Velocity  Increment  ft/aec 

33519 

29838 

Energy 

106  ft2/sec2 

561,8 

445.2 

Initial  State 

Radius 

I06ft 

20,93 

20.93 

Velocity 

ft/aec 

1492 

1492 

Potential  Energy 

10®  ft2/eec2 

-672.7 

-672.7 

Kinetic  Energy 

10‘ 

1.1 

1.1 

Total 

10°  ftVsec2 

•671.6 

-671.6 

-  il 


Final  State 


&  ' 

• 


Radius 

I06  ft 

30.18 

23.11 

Velocity 

Potential  Energy 

it/***  * 

10?  ftvsec* 

2J598 
-466  .4 

24682 

-609.2 

5 

Kinetic  Energy 

10®  ft2/sec2 

233.2 

304.6 

Total 

10®  ft*/ see2 

-233.2 

-304.6 

t 

Energy  Loss 

Jfif 

*  * 

Absolute 

10®  ft2 /sec2 

123.4 

78.2 

4 

Percent  of  Input 

22.0 

17.6 

In  order  to  show  the  effect*  of  specific  impulse  and 

* 

m*M  fraction,  it  teess  reoaonabla  to  select  a  total  ideal 
rocket  velocity  increment  which  ie  represents  five  of  that 
required  for  circular  orbit  entry  using  the  16  inch  gun  et 
Barbados.  This  has  been  done  for  a  velocity  of  24,400  ft/sec 
and  the ‘result*  are  shown  in  fig.  14.  Geometric  staging  was 
assumed. 

* 

Typical  trajectory  characteristics  are  shown  in  Pig.  15  * 

(height  vs.  range).  Fig.  16  (height  vs,  airspeed)  end  Fig.  17 

.  ** 

(height  ar.d  airspeed  vs.  time).  Typical  vehicle  parameters 
are  Hated  below. 


STAGE 

Propellant 

Total 

l 

1152 

288 

1440 

2 

320 

80 

400 

3 

96 

24 

120 

Payload 

0 

40 

40 

3  &  Payload 

96 

64 

160 

2,  3  &  Payload 

416 

144 

560 

l.  2,  3  &  Payload 

1568 

432 

2000 

12 


C0WCLP8I0KS  AMD  KEC^tgNPATIONS 

The  orbital  capabilities  of  a  multi-stage  gun  launched 
rocket  have  been  demonstrated.  In  particular,  it  has  been 
shown  that  with  the  existing  16  inch  gun  installation  at 
Barbados,  launching  a  2000  lb.  three  stage  vehicle  at  4500 
ft/sec.,  a  payload  of  40  lb.  can  be  placed  in  a  400  ran 
circular  orbit  (based  on  a  specific  impulse  in  vacuum  of  300 
sec.  and  a  propellant  mass  fraction  cf  0,8) 

Optimum  trajectories  fer  entry  into  circular  orbits  are 
close  to  the  g-t urn  type,  with  gun  elevation  in  the  neighbor¬ 
hood  of  45°,  first  stage  ignition  in  the  25-100,000  it.  range 
and  second  stage  ignition  as  soon  as  possible  after  first 
stage  burnout.  Optimum  vehicles  involve  third  stage  velocity 
Increments  somewhat  less  than  one- third  of  the  total  (for 
three  stage  vehicles)  hut  the  penalty  of  equal  velocity 
Increment  staging  is  small. 

Further  work  in  this  ares  could  be  directed  toward 
heavier  payloads  (:r  lower  specific  impulse,  or  lower  mass 
fraction),  lower  muzzle  velocity,  long  burning  times,  higher 
gun  elevations,  and  variations  in  drag  (or  atmospheric 
properties),  launch  location,  and  azimuth. 

It  is  hoped  that  the  optimization  procedure  will  be 
programmed  in  the  near  future,  so  that  t*-,e  optimum  vehicle 
and/or  trajectory  will  be  produced  automat ically,  at.  least  for 
s  limited  number  ->f  parameters. 


■w  «  '  *• 
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F  I  CURE  2 


ATMOSPHERIC  PRESSURE  APPROXIMATION  ERROR 
Bated  on  U.S.  1962  Standard  Afroaphere  (Ref,  5) 


FIGURE  3 
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SONIC  SPEED  APPROXIMATION  ERROR 
Based  rn  t’.S.  1962  Standard  Atmosphere  (Re£„  5) 
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FIGURE _ 5. 


EFFECT  OF  SECOND  STAGE  IGNITION  POINT 
Payload  16  Lb. 

Third  Stage  Weight  64  Lb.  (Geometric  Staging) 
First  Stage  Ignition  at  100,000  Ft. 

First  &  Second  Stage  Thrust  Direction  Type  0 


Elevation 


FIGURE 


EFFECT  OF  FIRST  STAGE  IGNITION  POINT 

Pay  I  cad  16  Lb 

Third  Stage  Weight  64  Lb.  (Geometric  Staging) 
Second  Stage  Thrust  Direction  Type  0 

Circular  Ocbit  Height 


First  Stage  TcriSt  Dlrec tic." 

-  Type  1 

•  Type  0 


Second  Stage  Igritioo.  Delay 
_____  0  Sec. 

X  10  Sec. 


Ga->.  Elevati 
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FIGURE  8 


EFFECT  OF  SECOND  STAGE  IGNITION  DELAY 


AND  fflRUST  DIRECTION 


Payload  40  Lb. 

Third  Stage  Weight  140  Lb. 

First  Stage  Ignition  at  25000  Ft. 
First  Stage  Thrust  Direction  Type  1 
Second  Stage  Thrust  Direction  Type  2 
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Path  at  Ignition 
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FIGURE  10 


EFFECT  OF  THIRD  STAGE  WEIGHT 

Payload  40  Lb. 

First  Stage  Ignition  at  25,000  £t. 
Second  Stage  Ignition  Delay  10  Sec. 
First  Stage  Thrust  Direction  Type  1 
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EFFECT  OF  OPTIMIZING  THIRD  STAGE  WEIGHT 


First  Stage  Ignition  at  25000  Ft. 
Second  Stage  Ignition  Delay  10  Sec. 
Firsr  Stage  Thrust  Directior  Type  1 
Second  Stage  Thrust  Direction 
Type  0  for  Payloads  <40  Lb. 
Optimized  for  Payloads  >40  Lb,, 
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FIGURE  12 


CQMTARIS1  *N  OF  NOMINAL  AND  OPTIMIZED 
THIRD  STAGE  WEIGHT 


First  Stage  Igntfl-i  at  25000  F*. 

Set  Stage  Delav  10  Sec 

First  Stage  Thrust  DirecM)'*  Typ*  1 
Sec  nd  Stage  TVust  l-lrectir'i  Tvre  0 


00 


50 


00 
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FIGURE  13 


EFFECT  OF  SECOND  STAGE  THRUST  DIRECTION 


Payload  50  Lb. 

Third  Stage  Weight  130  Lb. 

First  Stage  Ignition  at  25000  Ft. 

Second  Stage  Ignition  Delay  10  Se<-. 

First  Stage  Thrust.  Direction  Type  1 
Absolute  Flight  Pa th  Angle  at  Second  Stage 
Ignition  -  24  Deg. 


Gun  Elevation 
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FIGURE  14 


EFFECT  OF  SPECIFIC  IMPULSE  AND  MASS  FRACTION 


Rccket  Total  Ideal  Velocity  Ir.creTient 


24,400  ft/sec 


FIGURE  15 


TYPICAL  TRAJECTORY  -  HEIGHT  VS.  RANGE 

Payload  40  Lb . 

Specific  Impulse  300  See. 

Haas  Fraction  0.8 
Burning  Tinea  20  Sec. 

Gun  Elevation  39  Deg. 

Firat  Stage  Ignition  at  25,000  Ft. 

Firat  Stage  Direction  Type  1 
Second  Stage  Ignition  Delay  10  Sec. 

Second  Stage  Direction  Type  0 


Circular 


RANGE  INCREMENT 
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APPENDIX  I 


LIST  OF  SYMBOLS 


FORTRAN 

MATH 

DEFINITION 

AP 

P 

Atmospheric  Pressure 

AT 

0 

Latitude  (See  Note  (2)  ) 

ATD 

0° 

Latitude  (Degrees) 

ATG 

•1 

Latitude  of  Gun  (Degrees) 

ATO 

*8 

Latitude  of  Gun 

AVU 

Length  of  Unit  Vector  on  U  Axis 

Minus  1. 

AVW 

Length  of  Unit  Vector  on  W  Axis 

Minus  1, 

AX 

a 

Semi-Major  Axis  of  Orbit 

AZ 

fi 

Azimuth  (See  Note  (1 )  ) 

AZD 

e# 

Azimuth  (Degrees) 

AZG 

Azimuth  of  Gun  (Degrees) 

AZS  (N) 

fit 

Azimuth  of  Thrust  Vector  at  Ignition 
(Degrees) 

C  (I) 

h 

Kutta-illl  K  for  Variable  I 

CAZ 

C08 

fi 

CD 

cd 

Drag  Coefficient 

CEL 

COS 

s 

CLA 

COS 

0 

CNA 

Normal  Force  Coefficient-Angle  of 
Attack  Derivative 

CT 

c 

Iteration  Count 

DL 

Computing  Time  in  60tbs  of  a  Second 

DLS 

Computing  Time  in  Seconds 

DR 

D 

Drag 

•  »  OH’gWSSfSf 


I 


••C-WnjSi 


»wi  1ml  ii 
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Li;  U 

MATH 

DEFINITION 

DT 

At 

Step  Size 

DTB  (N) 

> 

ft 

O* 

Burning  line 

DTD 

AtP 

Possible  Step  Size;  Previous  Step 

Size 

DTI 

Atx 

Noninsl  Step  Size 

DTP 

Atp 

Nominal  Print-Out  Interval 

DU 

U-Directioc  Cosine  of  Thrust  Vector 
(Body  Axis) 

DV 

«v 

V-Direction  Cosine  of  Thrust  Vector 
(Body  Axis) 

OH 

Sw 

W-Directlon  Cosine  of  Thrust  Vector 
(Body  Axis) 

OX 

s* 

X-Dlrection  Cosine  cf  Thrust  Vector 
(Body  Axis) 

DY 

Sy 

Y-Direction  Cosine  of  Thrust  Vector 
(body  Axis) 

DZ 

% 

Z-Direction  Cosine  of  Thrust  Vector 
(Body  Axis) 

EA  (N) 

A 

Exhaust  Area 

KL 

• 

Elevation 

KID 

•°  . 

Elevation  (Degrees) 

ELG 

i° 

•i 

Elevation  of  Gun  (Degrees) 

ELS  (N) 

iS 

Elevation  of  Thrust  Vector  at  Ignition 
(Degrees) 

m 

M 

Mach  Number 

r 

P 

108/TFf 

PN 

Normal  Force 

FNU 

*nu 

U-Coaq>ocant  of  Normal  Force 

FNV 

^nv 

V- Component  of  Nonna l  Force 

FEW 

^nw 

W-Component  of  Nv/rmal  Force 

"v'r ' *-*ry 4*-* 
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FORTRAN 

MATH 

DEFINITION 

FR  (N) 

f 

Mass  Fraction 

GC 

ke2 

Gravitational  Constant 

GO 

8o 

Acceleration  due  to  Gravity  at  S.L. 

HA 

ha 

Apogee  Height 

HP 

«P 

Perigee  Height 

HT 

H 

Height  (Feet) 

HTN 

Height  (Nautical  Miles) 

I 

1 

Variable  Number 

11  (N) 

h 

Ignition  Indicator 

12  (N) 

l2 

Firing  Indicator  (Direction) 

J 

N 

Integration  Step 

K 

Error  Indicator 

L 

Computer  Clock  Reading  at  End  of  Run 

LO 

Computer  Clock  Reading  at  Beginning 

of  Run 

M 

Case  Number 

N 

n 

Stage  Number 

NL 

nL 

Number  of  Lines  Printed  this  Page 

NN 

nl 

Previous  Stage  Number 

NP 

n 

P 

Page  Number 

NS 

n8 

Number  of  Stages 

0N 

e 

Longitude  (See  Note  (4)  ) 

0ND 

8° 

Longitude  (Degrees) 

0NG 

•S 

Longitude  of  Gun  (Degrees) 

0NO 

®g 

Longitude  of  Gun 

PA 

1 

Absolute  Flight  Path  Angle  (See  Note 

(2)  ) 

i 


i 


*&A 
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FORTRAN 

MATH 

DEFINITION 

FAD 

V • 

Absolute  Flight  Psth  Angle  (Degrees) 

FAF 

y* 

Absolute  Flight  Path  Angle  (Previous 
Step) 

Q  (I) 

Qi 

Kutta-Gill  Q  for  Variable  I 

QN 

QS 

QS 

QI  00 

On 

Ignition  Value  (See  Note  (5)  ) 

RA 

d 

Range  (Feet) 

RAN 

Range  (Nautical  Miles) 

RE 

re 

Earth  Radius 

RF 

rP 

Perigee  Radius 

SA  (N) 

S 

Cross  Section  Area 

SAZ 

sin  0 

SEL 

•In  ft 

SI  00 

*«0 

Specific  Impulse  in  Vacuum 

SL 

P 

Semi • La t us  Rectum  of  Orbit 

SLA 

sin  0 

SF 

Jl 

Angular  Rate  of  Earth 

sywu 

Scalar  Product  of  Unit  Vectors  on  U 
and  W  Axes 

ss 

c 

Speed  of  Sound 

S3 

s3 

Event  Indicator 

S4 

S4 

Firing  Indicator  (Condition) 

T 

t 

Time 

TB 

fB 

Time  of  Burnout 

TH 

T 

Thrust 

THO 

T- 

Thrust  in  Vacuum 
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FORTRAN 

NATH 

DEFINITION 

TI 

*1 

Time  of  Ignition 

T# 

lo 

Time  rf  Last  Print-Out 

V 

V 

Absolute  Velocity 

VR 

vr 

Relative  Velocity  (Referred  to  Earth) 

VRU 

Vru 

U-Component  of  Relative  Velocity 

VRV 

Vrv 

V-Component  of  Relative  Velocity 

VRW 

Vrw 

W-C mp-'nent  of  Relative  Velocity 

VX 

vx 

X-C omponent  of  Unit  Vector  on  V  Axis 

VY 

vy 

Y-fomponent  of  Unit  Vector  on  V  Axis 

VZ 

vz 

Z-Cjmponent  of  Unit  Vector  on  V  Axis 

WAZD 

S* 

Azimuth  of  Relative  Velocity 

WELD 

®w 

Elevation  of  Relative  Velocity 

WI  (N) 

W 

Total  Weight  Before  Ignition 

WU 

Wu 

U-Component  of  Angular  Velocity  Vector 

Y  (I) 

Yi 

Variable  I 

YD  (I) 

Yi 

Time  Derivative  of  Variable  I 

NOTES: 

1.  Azimuth  *8  measured  from  the  gun  (the  input  gun 

direction)  is  measured  in  a  horizontal  plane, 
clockwise  from  North  as  viewed  from  above  (farther 
from  the  earth  center). 

Azimuth  as  measured  f rom  the  vehic ie  (the  input 
direction  of  the  thrust  vector  at  ignition  and  the 
output  bcdy  axis  and  wind  directions)  is  measured 
in  a  horizontal  plane,  clockwise  from  the  horizontal 
absolute  velocity  vector  as  viewed  from  above. 


2. 


Elevation  and  flight  path  angle  are  measured  from  the 
horizontal  plane,  positive  upwards. 
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MOTES  (Cont’d) 

3.  Latitude  la  measured  positive  North  from  the 
equator. 

4.  longitude  la  meaaured  positive  Eaat  from 
Greenwich. 

5.  The  ignition  value  la  interpreted  aa 

Height  In  thouaanda  of  feet  if  Ij  ■  0 

Flight  path  angle  in  degrees  if  1^  ■  1 

lime  after  previous  burnout  if  1^  ■  2 


BASIC  VARIABLES 


VARIABLE  NO. 

I 

MATH 

SYMBOL 

DEFINITION 

l 

W 

Weight 

2 

cf 

Radius  Rate 

3 

r 

Radius 

4 

<*w 

W-Component 

of 

Angular  Velocity 

5 

ux 

X-Component 

of 

U-Direction 

6 

Uy 

Y-Component 

of 

U-Direction 

7 

Z- Component 

of 

U-Direction 

8 

X- Component 

of 

W-Direction 

9 

Wy 

Y-Component 

of 

W-Direction 

10 

w 

"a 

Z-Component 

of 

W-Direction 

Variable  Y  (I) 

Rate  of  Change  of  Variable  YD  (I) 


Mm***  ♦*->  -K^-rt^-*****^ 
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LIST  OF  SYMBOLS  (Cart'd) 
INDICATORS 


INDICATOR _ VALUE _ MEANING _ 

xl  0  Ignition  at  Preset  Height 

1  Ignition  at  Preset.  Flight  Path  Angle 

2  Ignition  at  Preset  Time  After  Burnout* 
NOT  0,1,2  Igritior  at  Burnout* 

12  NOT  1,2,3  Firing  at  Constant  Attitude  Initially 

larger*  *o  Trajectory 

1  Firing  in  Wind  Direction 

2  Firing  at  Preset  Attitude 

3  Firing  Along  Trajectory 

S3  -1.  Impact 

0  Burn: :  t 

1.  IgMtlon  or  Launch  or  Apogee 

2.  No  Evert 


S4  -1.  No*.  Firing 

0  No  Mere  Stages  to  Fire 

1.  Firing 


*  Burnout  refers  to  the  previous  stage,  or  in  the  case  of  first 
xtage  ignition,  the  muzzle. 


L-, -I  "5  t  5*-*  ' 


APPENDIX  II 


MATHEMATICS 
Initialise  Run 


VR  sin  S 
re  +  H 

£  (VR  co*  8  sin  £  +«/Xr  co*  0)^  +  (VR  co*  8  co*  {5  )2J  ^ ft 

CO*  0 


0 

•in  0 


-  (VR  cos  f  sin  6  +  JT  r  cos  0)  sin  0/ru£ 

-  VR  cos  E  cos  fi  I*  OK, 

MK  W 

(VR  cos  •  sin  <3  +_fl  r  coj  0)  cos  #/rfcrw 


Initialise 


Vx  - 

UpWy  -  UyWg 

V  - 

u  w  -  u  w 

y 

x  •  ax 

v.  - 

U„Ww  -  u  w 

• 

y  x  x  y 

VRU  * 

4 

■ 

> 

r (Ww  -ilV,) 

VEW- 

rJl  Vg 

vR  - 

(vt02  +  V,v2  +  v^2 

H  - 

r  -  re 

V  - 

C4  2  +  <ru^)2J  * 

V  » 

tan*1  (^/rwj,) 

At1- 

(q/f  .|)At  /<*-* 

(1000Q  -  H)/*( 

■ 

Q  -  t  +  tB 

■ 

-  H/<( 

K 

-d/i 

t  ^  tB  n*ng 


n<ng  It  -  2 
n  *  ng  ^.0 

n  *  ng 


~4 


-f' 


,  i- v-», ,  vy(?m*jltv*.-aa>  vr '  1-  ? 


2  - 


®u 

m 

*/v 

) 

) 

) 

) 

) 

®v 

m 

r*w/V 

x2 

H  1,  2  or  3 

K 

m 

0 

*u 

m 

sin  8 

) 

) 

•v 

m 

cos  8  cos  P 

) 

) 

l2 

-  2 

*w 

m 

-  cos  8  sin  (3 

) 

sx 

m 

8UUX  +  fivVx  +  8wWx 

) 

) 

) 

) 

®y 

m 

8uUy  +  «vVy  ♦  ®„Wy 

*2 

■  1  or  3 

*« 

m 

*uug  ♦  svv«  +  SwW* 

) 

• 

w 

m 

f  <»n  +  1  "  V  fcB 

To 

m 

-  IW 

0 

m 

tan-1  (Uy/Ux)  -Jit  +  0O 

0 

m 

“"'‘["A’  +  Uy2>  *] 

Calculate  Variables 


V* 

• 

U.Wy  -UyW# 

vy 

m 

Ux«a  -  U*Wx 

v. 

m 

uywx  -  uxwy 

m. 

m 

* 

VRV 

m 

r  (urw  -JlVm) 

VRW 

m 

rJLVa 

VR 

m 

(y  ^+U  ^  +  V 

H 

m 

r  ’  re 

P 

m 

2116. 22e  ’4‘42  x  10‘  H 

H  106 

P 

m 

0 

h  y  106 

VtWM&f'lk 


i 


!: 

I: 


i 

\ 
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Calculate  Thrift 


T 

m 

T0  -  AF 

•u 

m 

«*Ux  *  W  ♦  «aUa 

4 

®v 

m 

Vx  ♦  Vy  +  ®»v« 

• 

I2  U  or  3 

K 

m 

•xwx  +  V'y  +  *«wa 

a 

•u 

m 

v*u/vi 

*v 

m 

vt»At 

* 

*2-  l 

r% 

Vw/vt 

•u 

m 

f /v 

* 

®v 

m 

► 

1 

X2“3 

•* 

m 

0 

4 

Calculate  Drat  and  Normal  Force 

N  -  Vt/c 

Cd  -  0.35  M<1 

Cd  -  0.18  ♦  a’M/2  M^l 

QB  •  0.7  m2  S 

D  *  Cd  Q8 

-h/4  „m 

rn  2  +  12  e  -  10  e 

On  *  cm  Q8/Va 

*nu  *  tC*u  ®vv»v  +  *»VW>  *  vtu  <«v2  +  »V2)J 

fw  ’  In  C,  <•«  *»u  +  V»)  -  v*v  <«u2  ♦  s«23 

fnw  -  %  C*»  <»u  VRU  +  «,Vrv>  -  vw  «u2  +  *v2j] 

rn  *  »nu2  + 'nv2  + 'w2>* 


•'it "" 
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Calculate  Rates 


sr 

- 

rurw2  +  (g0/M)  (Tfiu  - 

• 

r 

- 

«"w 

- 

[  -2  {Ufc  +  (  g  0/W)  (T8, 

- 

(ISU  -  D  Vm/VR  +  Fnw) 

Ux 

SI 

VxWw 

Uy 

r 

U. 

X 

< 

N 

Wx 

- 

-  V  Ci/' 

X  u 

Wy 

- 

*  ^yt*,u 

W. 

* 

-  vzK. 

Output 


-1  <uy/ux) 

-ilt  +  e0 

_l  <“./  1/ 

»X2  +  »y2> 

i*1  [sin  (8 

-  eo)  cos  (0  -  0o)/sin  (0 

-  #00 

tan-1  Qsin  (8  -  80)  1  / 1  sin  cos 

(0  -  », 

"l  <VRU  / 

VVRV2  +  VRH2> 

'*  (  -  Vrw 

/  VRV) 

tan-1  <SU  /  Vsv2  +  8„2) 
tan'1  (-  «w/*v) 
l/ (2/r  -  V2/ke2) 

(Vrcos  V)2  /ke2 
p  /(l  +  1-p/a) 


rP  '  re 


)  Not  firing 


)  Firing 


ap/rp  -  re 


Inf  cate 

Let  X  represent 
For  ell  Steps 

Step  1 

Step  2 

Step  3 

Step  4 

After  each  step, 


5  - 


W  ^  r  U rv  ux  Uy  U#  wx  Wy  or 

Ax  ■  x  A tj 

x*  •  x  +  Ax/2 
A  -  Ax 
t  -  t  +  A  tx/2 

x1  -  x  +  (l  -  1/  Jl)  (Ax  -  A  ) 

A  -  (2  -  Si)  Ax  +  (  -2  +  3/  «T2)  A 

X1  -  X  +  (1  +  1/  Si)  (Ax  -  A) 

A  -  (2  +  Si)  Ax  +  (  -  2  -  3/lT”2)A 
t  -  t  +  Atj/2 

*l  -  x  +  A  x  /6  -  A  /3 


x  is  recalculated  for  the  next  step. 


APPENDIX  III 


COMPUTER  PROGRAM 


(SOURCE  DECK  LISTING) 

FORTRAN  IV  -  IBM  7040 
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H**P  TRAJ.  PROG.  R.  M.  MCKEE  FORTRAN  SOURCE  LIST  R MCKEE 

1SN  SOURCE  STATEMENT 

C 

C  HARP  TRAJECTORY  PROGRAM  -  R.  M.  MCKEE 

C 

C  INPUT 

'  c 

1  DIMENSION  AZS( 10 ) ,C( 10 ) ,DTB( 10 ) ,ELS(1Q) ,£A( 10) ,FR( 10) , 

1  11(10),I2(10),Q(10),QT(10),SA(10),SI(*9),WT(10)«Y(10),YD(10) 

2  1000  CAU  CLOCK  (L) 

3  RE AO (St 1001 )M,NS, ( WT (NS+1 ) ) ,  ( SAINS+l ) )  ,OTI ,DTP, 

1  VRtELGt AZGtONG* ATGtHT 

16  1001  FORMAT  (2I5,F5.2,F5.3,3F5.0,4F5*2,F5.0) 

1?  IF  (NS.EQ.O)  GO  TO  2 

22  READ  (5*1002)  ( I 1(N> * I2(N) ,WT(N) tSA(N) «EA(N) ,DTB(N) * OT (N) ♦ 

1  ELS(N),AZS(N),FR(N»tSI(N),N«l,NS) 

27  1002  FORMAT  (215* F5.lt2F5.3t 5F5.2* F5.0) 

30  2  NP  >  1 

31  WRITE  (6*1003)  M*NP 

32  1003  FORMAT  ( 1HI, 15X.41HHARP  TRAJECTORY  PROGRAM  *  R.  M.  MCKEE* 

1  35X,4HCASE,I5,5X,4HPAGE,I3) 

33  WRITE  (6* 1004 )M* VR 

34  1004  FORMAT  I1H-.10X,  UHCASE  NUMBER. 6X. I7.9X, 15HNUZZLE  VELOCITY* 

1  2X.F7.0) 

35  WRITE  (6, 1005 )NS*ELG 

36  1005  FORMAT  ( 1H  .10X.13HNC.  OF  STAGES, 4X, I7.9X,  9HELEVATI0N*6X«F7.2) 

37  WRITE  (6* 1006 )WT (NS+1 ) , AZG 

40  1006  FORMAT  < 1H  * 10X* 15HPAYL0A0  WEIGHT, 2X.F7. 2, 9X,  7HAZIMUTH*10X,F7.2) 

41  WRITE  (6, 1007 ) SAINS'*' 1 )  »ONG 

42  1007  FORMAT  ( 1H  , 19X,4HAREA,4X,F7.3.9X,9HL0NGITU0E,8X,F7.2) 

43  WRITE  (6.1008)OT!*ATG 

44  1008  FORMAT  ( 1H  ,10X,  9HSTEP  SIZE*8X,F7.0,9X*8HLATITU0E,9X*F7.2) 

45  WRITE  (6* 1009)DTP»HT 

46  1009  FORMAT  ( 1H  , 10X* 15H0UTPUT  INTERVAL, 2X,F7.0,9X,6HHEIGHT,11X,F7.0) 

47  IF  (NS.EQ.O)  GO  TO  1 

52  WRITE  (6,1010) (N.N»1, NS) 

57  1010  FORMAT  ( IH-* 10X * 12HSTAGE  NUMBER, 6X, 6115) 

60  WRITE  (6,10il)(Il(N),N«l,NS) 

65  1011  FORMAT  (1H0,10X*18H IGNITION  INDICATOR, 61 15 ) 

66  WRITE  (6,1012)(I2(N),N«1,NS) 

73  1012  FORMAT  ( IH  , 10X , 16HF IR ING  IN0ICAT0R.2X.6I 15 ) 

74  WRITE  (6»1013)(WT(N),N*1,NS) 

101  1013  FORMAT  ( IH  * 10X, 14HWEIGHT  ( TOTAL ) ,4X,6F15. 1 ) 

102  WRITE  (6,l014)(SA(N),N«l,NS) 

107  1014  FORMAT  ( IH  , 10X .4HAREA , 14X ,6F 15. 3 ) 

110  WRITE  (6,1015) (EA(N),N«1, NS) 

115  1015  FORMAT  ( IH  , 10X,  12HEXHAUST  AREA,6X,6F15.3) 

116  WRITE  (6,1019) (OTB(N),N«l, NS) 

123  1019  FORMAT  ( IH  , 10X , 12HBURN1NG  T 1  ME , 6X , 6F 1 5. 2 ) 

124  WRITE  (6*i016)(QT(N)*N»l,NS) 

131  1016  FORMAT  ( IH  * 10X, 14HIGNITIGN  VALUE ,4X,6F15.2 > 

132  WRITE  (6,1017)(ELS(N),N«1,NS) 

137  1017  FORMAT  (IH  , 10X, l 6HF IR ING  ELEVATION, 2X.6F15.2) 

140  WRITE  (6,1018)  (AZS(N),N»l, NS) 

145  1018  FORMAT  ( IH  , 10X, 14HF IR ING  AZIMUTH, 4X,6F15. 2) 

146  WRITE  (6,1026)  ( FR( N ) , N»l , NS ) 

153  1026  FORMAT  ( IH  ,10X,13HMASS  FRACT ION, 5X,6F 15.2 ) 
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15*  WRITE  16» 1020) (SI  (N),N»l,NS) 

161  1020  FORMAT  ( 1H  * 10X, 16HSPECIFIC  IMPULSE, 2X,6F15.0) 

162  1  WRITE  (6,1025) 

163  1025  FORMAT  (1H-,  7X,*HTIME.6X,6HHE1GHT,7X,5HRANGE,*X,8HVEL0CITY, 

1  *X,8HAIRSPEE0,2X,10HPATH  ANGLE, 3X,9HELEVATI0N»5Xt 7HAZIMUTH, 

2  5X,7HWIND  EL,5X,7HWIND  AZ/1H  ) 

16*  NL  «  28 

165  IF  (NS.EQ.O)  NL  *  13 

C 

C  INITIALIZE  RUN 

C 

170  10  GC  ■  1.A07639E16 

171  GO  «  32.1*65 

172  RE  «  2.09256AE7 

173  SP  •  7.31958  E  -5 

17*  F  *  57.295780 

175  ONO  «  ONG/F 

176  ATO  «  ATG/F 

177  AZ  »  AZG/F 

200  EL  «  ELG/F 

201  N  «  1 

202  CT  «  0 

203  S*  »  -1. 

20*  IF  (NS.EQ.O)  S*  «  0 

207  UT  -  DTI 

210  DTD  *  0 

211  TH  «  0 

212  YO ( 1 )  *  0 

213  T  *  0 

21*  TO  *  0 

215  TI  *  0 

216  TB  *  0 

217  Y 1 1 )  *  WT  1 1 ) 

220  SLA  »  SiN(ATO) 

221  CLA  «  COS ( ATO ) 

222  SAZ  «  SIN(AZ) 

223  CAZ  *  COS ( AZ ) 

22*  SEL  -  S ! N ( EL ) 

225  CEL  *  COS ( EL ) 

226  Y ( 2 )  =  VR*SEL 

227  Y ( 3 )  »  RE+KT 

230  Y (  * )  *  SQRT((VR«*CEL*SAZ*SP*Y(  3 ) »CLA ) **2+ ( VR«CEL«CAZ )*»2 )  /  V(3> 

231  Y ( 5 )  *  CLA 

232  Y ( 6 )  *  0 

233  Y(7)«  SLA 

23*  Y ( 8 )  *  -(VR*CEL*SAZ+SP*Y(3)«CLA)«SLA/(Y(3)*Y(A)) 

235  Y ( 9 )  «  -VR*CEL*CAZ/(Y(3)*Y(*)  ) 

236  Y ( 10 )  *  ( VR*CEL*SAZ+SP*Y ( 3 ) *CLA ) *CLA/ l Y ( 3 ) *Y(* I ) 

C 

C  INITIALIZE  STEP 

C 

237  20  VX  *  Y(7)*Y(9)-Y(6)*Y( 10) 

2*0  VY  *  Y(5)*Y(10)-Y(7)«Y(8) 

2*1  VZ  *  Y(6)*Y(8)-Y(5)*Y(9) 

2*2  VRU  *  Y ( 2 ) 
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A 


243 

244 

245 

246 

247 
250 

253 

254 

255 
260 

263 

264 

265 
270 
273 
276 
301 
304 
307 
312 
315 

320 

321 

324 

325 

326 

327 

332 

333 

334 

335 

336 
341 

344 

345 

346 
34  7 

350 

351 

352 

353 

354 

355 

356 

357 

360 

361 

362 

363 

364 

365 

366 

367 
370 


VRV  «  Y( 3 ) • ( Y( 4 )-SP*Y ( 101) 

VRW  «  Y( 3 )«SP*VZ 

VR  ■  SQRT(VRU»*2+VRV«*2*VRW*«2) 

HT  »  Y (31-RE 

V  -  SORT (Y(2)»*2'MY(3J»Yt4)  )*«2 ) 


IF 

(T.NE.TB)  PAP 

M 

PA 

PA 

«  ATAN2 ( Y ( 2 ) , 

Y  ( 3 

) *Y  ( 4 ) ) 

S3 

«  2. 

IF 

(T.EQ.O.)  S3 

«  l 

• 

IF(HT.LT.0.)S4»0 

• 

OTO 

1  *  OT 

OT 

■  OTI 

IF 

( S4.EQ.0. )G0 

TO 

27 

IF 

(S4.E0.1.)  GO 

TO 

24 

IF 

(Il(N).EQ.l.ANO. 

T.EQ. 

TB) 

GO 

TO 

30 

IF 

(Il(N).EQ.l) 

OTO  • 

( QT(N) 

/F-PA)*DTD/( PA-PAP) 

IF 

(Il(N).EQ.O) 

OTD  « 

(1000. 

•QT(N)-HT)/Y(2) 

IF 

( I1(N).EQ.2 

) 

DTD  ■ 

OT  (Nl- 

T+TB 

IF 

( I 1 (N) .EQ.3) 

GO 

TO  21 

IF 

(ABS(OTO).GE. 

OTI 

)  GO 

TO 

30 

IF 

(ABS(OTO).LT. 

.01 

)  GO 

TO 

22 

CT 

-  CT+1. 

IF 

(CT.GT.20.) 

GC 

1  TO  29 

OT 

«  OTO 

GO 

TO  30 

22  CT 

«  0 

IF 

(T.GE.TB)  GO 

TO 

21 

NN 

■  N-I 

WRITE  (6,1107)  N,NN 

1107  FORMAT  ( 1H0, 24HATTEMPTED  IGNITION  STAGE, 12, IX, 
1  20HBEF0RE  BURNOUT  STAGE, 12) 

GO  TO  70 

21  IF  ( I2(N).EQ.l.0R.I2(N).£Q.3)  GO  TO  23 
IF  (I2(N).EQ.2  )  GO  TO  25 

OU  «  Y (21/V 
OV  »  Y ( 3 ) *Y ( 4 ) /V 
OW  »  0 
GO  TO  26 

25  EL  »  ELS(N)/F 
AZ  «  AZS(N)/F 
CEL  »  COS( EL ) 

SEL  «  SIN ( EL ) 

CAZ  «  COS(AZ) 

SAZ  -  SIN(AZ) 

OU  »  SEL 

OV  «  CEL*CAZ 
OW  «  -CEL*SAZ 

26  OX  «  DU«Y (5)+DV*VX*DW»Y(8) 

OY  »  OU«Y ( 6  >  +OV*VY+OW»Y ( 9 ) 

OZ  »  DU*Y ( 7 ) ♦OV»VZ*OW*Y ( 10 ) 

23  S4  *  l. 

S3  ■  1. 

YDt l )  «-FR(N)MWT(N)-WT(N*i) )/DTB(N) 

THO  —  SI(N)«Y0(1) 

T I  »  T 


*W«*w.*  ,, 


'  *•  tv»-  ‘•'W!*  ' 
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1101 


1103 


GO  TO  30 


OT  «  DTB(N)/4. 

IF  (DT.GT.DTII  OT  *  OT I 

0ND«  F*(ATAN2(Y(6),Y(5))-SP«T«-0N0) 

ATD«  F»(ATAN2(Y(7)» SORT (Y(5)**2+Y(6)»*2)  )  ) 

WRITE  (6,  HOl'/NtONOf  ATO 

FORMAT  ( 1H  . 5HST AGE ,12, IX t20H IGNITION  -  LONGITUDE, F9. 3* 
1  9H  LATITU0E,F9.3) 

NL  *  NL  *1 
GO  TO  30 

IF  (  (T+OT).GE.(TI+DTB(N>))  OT  «  TI*DTB(N)-T 
IF  (0T.LT..01)  S3  »  0 
GO  TO  30 

IF  ( Y( 2 ) «GE«0« )  GO  TO  102 
OTO  «  -HT/YI2) 

IF  (ABS(OTO).GT.DTI)  GO  TO  30 
IF  ( ABS(OTD) . LT. .01 )  GO  TO  28 
CT  *  CT  +  1. 

IF  (CT.GT.20.)  GO  TO  29 
OT  «  OTO 
GO  TO  30 

ONO*  F* ( ATAN2 ( Y ( 6 ) ,  Y(5) )-SP*T+ONO) 

ATO*  F»(ATAN2(Y(7)»  SORT  CYC  5) *»2^Y 16) **2 ) I ) 

WRITE  (6,1103)  ONO, ATO 

FORMAT  (1H  , 18HIMPACT  -  LONGITUDE, F9. 3, 9H  LATITUDE, F9. 3) 
NL  *  NL+1 
S3  «  -l. 

GO  TO  30 

IF  (T.EQ.T8)  GO  TO  30 
OTO  *  -Y(2)/Y0(2) 

IF  ( ABS (OTO). GT .DTI )  GO  TO  30 
IF  (ABS(0T0).LT..01)  GO  TO  103 
CT  »  CT+1. 

IF  (CT.GT.20.)  GO  TO  29 
OT  «  DTD 
GO  TO  30 


1109 


29 

1104 


WRITE  (6,1109) 

FORMAT  (1H  ,6HAP0GEE ) 

NL  *  NL  +  1 
S3  »  1. 

GO  TO  30 
WRITE  (6,1104) 

FORMAT  (1H0,17HITERATI0N  FAILURE) 
GO  TO  70 


CALCULATE  VARIABLES 


30  00  65  J  «  1,4 

IF  (NL.LT.55)  GO  TO  31 
NP  *  NP+1 

WRITE  (6,1003)  M , NP 
WRITE  (6,1025) 

NL  »  5 

31  VX  *  Y(7)*Y(9)-Y(6)»Y( 10) 
VY  =  Y(5)*Y(iO)-Y(7)*Y(8) 


^  *-.t  i 
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504  VZ  »  Y(6)«Y(8)-Y(5)*Y(9) 

909  VRU  •  V ( 2 ) 

906  vrv  «  vm«m4)-sp*Y(ion 

907  VRW  •  Y( 3 )«SP*VZ 

910  VR  «  SQRT(VRU»»2*VRV*«2+VRW*«2) 

91!  HT  •  Y(3)-RE 

912  IF  (HT.LT.1E6)  GO  TO  32 

919  AF  >  0. 

916  OR  —Or 

917  GO  TO  40 

920  32  AP  >  2U6.22«EXP(-4.42E~5»HT) 

C 

C  CALCULATE  THRUST 

C 

921  40  IF  (S4.NE.  1.)  GO  TO  90 

924  TH  «  THO-EA(N)*AP 

929  IF  (I2(N).E0.1  )  GO  TO  41 

930  IF  (12(N).EQ.3)  GO  TO  43 

933  DU  •  DX*YI5)>*0V«V(6)«-0Z*V(7I 

934  OV  •  OX»VX  ♦OY*VY  +DZ*VZ 

935  OW  •  DX*V ( 8 )+DY»Y 19 1 ♦OZ»Y( 101 

936  GO  TO  90 

937  41  DU  «  VRU/VR 

540  OV  •  VRV/VR 

941  ON  *  VRU/VR 

942  GO  TO  50 

343  43  DU  •  Y12I/V 

544  DV  -  Y(3I*Y(4)/V 

945  OW  «  0 

C 

C  CALCULATE  DRAG  AND  NORMAL  FORCE 

C 

946  50  IF  (HT.LT. 36200.)  SS  •  1116.45  -  »0040986»HT 

951  IF  <MT. GE. 36200. .ANO. HT.LT. 65800. )  SS  •  968.08 

554  IF  (HT.GE. 65800. .ANO. HT.LT. 105000.)  SS  •  924.28  ♦  .0006658*HT 

557  IF  (HT.GE. 105000.. ANO. HT.LT. 155500. )5S  ■  811.54  ♦  .0017394*HT 

362  IF  (HT.GE. 155500. .ANO.HT.LT. 172000. )SS  •  1082.02 

565  IF  ( HT .GE. 172000. .ANO.HT .LT. 202000. ) SS  «  1291.40  -  .0012173*HT 

570  IF  (HT.GE. 202000.. ANO. HT.LT. 262500. ) SS  ■  1584.72-  .0026694«HT 

573  IF  (HT. GE. 262500. )  SS  »  884.0 

576  EM  «  VR/SS 

577  IF  (EM.LT.l.)  CO  «  .35 

602  IF  (EM.GE.l. )  CO  «  «18*EXP(~EM/2.) 

605  QS>.7«AP*EM*»2«SA(N) 

606  OR  >  CD*QS 

607  IF  (S4.NE.1..0R.I2(N).EQ.l)  GO  TO  51 

612  CNA  »  2.+12.»EXP(-EM/4.)-l0.»EXP(-EM> 

613  QN  •  CNA*QS/VR 

614  FNW  *  QN»(0U»(DV«VRV*0W»VRW)-VRU»(DV»*2^DN»»2) ) 

615  FNV  «  QN*(OV»(OU*VRU*OW*VRW)-VRVMDU«*2+OW«*2) ) 

616  FNW  »  QN*(DW»(DU*VRU*DV*VRV)-VRW*(DU»»2*DV**2) ) 

617  FN  ■  SORT (FNU»*2+FNV**2*FNW»«2) 

620  IF  FN/VI 15.LT.20.)  GO  TO  55 

623  WRITE  (6*1115) 

624  1115  FORMAT  (33H  NORMAL  ACCELERATION  EXCEEDS  20  G) 


FORTRAN  SOURCE  LIST  RMCKEE 


HARR  TRAJ.  FROG.  R.  M.  MCKEE 

ISN  SOURCE  STATEMENT 


625 

626 
627 
630 


51 


C 

C 

C 


C 

C 

C 


GO  TO  70 
FNU  »  0. 
FNV  »  0. 
FNW  «  0. 


CALCULATE  RATES 


631 

55  VDI2 ) 

Y 1 3 ) *V I 4 ) »»2* I GO/Y It) ) • I TH*DU-DR*VRU/VR*FNU )-GC/ VI 

632 

Y0(3) 

VI2) 

633 

Y0I4) 

I-2.»YI2)*YI4)^(G0/VI1))*(TH*0V-0R*VRV/VR*FNV))/V( 

634 

WU 

( TH*DW-DR* VRW/ VR +FNW )»GO/(Y(l)»Y(3)»V(4)) 

635 

YDI5) 

VX»V 14) 

636 

Y0C6) 

VY*V(4) 

637 

Y0(7) 

VZ« YI4) 

640 

YO  ( 8 ) 

-VX»WU 

641 

Y0I9) 

-VY»WU 

642 

YDI 10) 

•  -vz»wu 

OUTPUT 


643 

80 

IF  (S3,EQ.2..AN0.T.LT.!T0*DTPt.0R.J.NE.l)  GO  TO  60 

646 

TO  ■  T 

647 

ON  •  ATAN2(Y(6),Y(5))-SP*T*0N0 

650 

AT  «  ATAN2 I Y( 7) *SQRTI YI5)*»2*V16)*»2) ) 

651 

IF  (T.EQ.O.)  GO  TO  83 

654 

AZ  «  ATAN2I SIN(ON-ONO)*COSIAT-ATO) .SINIAT-ATO) ) 

655 

RA  «  RE*ATAN2(ABS(SIN(0N-0N0) ) »ABS(SIN(AZ) )»COSION-l 

656 

GO  TO  84 

657 

83 

RA  «  0 

660 

84 

RAN  >  RA/6076* 104 

661 

HTN  ■  HT/6076.104 

662 

HELD  «  F*ATAN2(VRU»SQRTIVRV»»2+VRW**2) ) 

663 

WAZD  «  F*ATAN2I-VRW.VRV) 

664 

IF  (S4.EQ.1.)  GO  TO  81 

667 

ELD  «  WELO 

670 

AZO  ■  WAZD 

671 

GO  TO  82 

672 

81 

ELO  «  F*ATAN2(0U.SQRTIDV*»2*DW«*2>) 

673 

AZO  «  F«ATAN2(-DW.0V) 

674 

82 

PAD  ■  F»PA 

675 

WRITE  16. l 105 ) T »HTN » RAN. V»VR» PAD. ELO. AZD. WELD. WAZD 

676 

1105 

FORMAT  (3F12.2.2F12.0.5F12.2) 

677 

NL  »  NL  +1 

700 

IF  (S3.NE.0.)  GO  TO  60 

703 

ONO  ■  F»ON 

704 

ATD  «  F*AT 

705 

WRITE  (6.  1102 )N.OND» ATO 

706 

1102 

FORMAT  (1H  . 5HSTAGE. I2»20H  BURNOUT  -  LONGITUDE. F9. 3 

1  9H  LATITUDE. F9.3) 

707 

NL  »  NL*l 

710 

N»N+l 

711 

S4  «  -l. 

712 

IF  (N.GT.NS)  S4  ■  0 

715 

YOU)  «  0 

716 

TH  «  0 

— — — ~ — .  J 


v 
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TIT  78  •  T 

T20  OT  •  DTI 

T21  Ylll  >  WT  ( N ) 

T22  IF  (HT.LT.2.5E5)  GO  TO  20 

723  AX  »  l./(2./Y(3)-V*«2/GC) 

T26  SL  •  (V«Y(3)*COS(PA))*«2/GC 

T2T  IF  (SL/AX.LT.l.)  GO  TO  85 

732  WRITE  (6,1110) 

733  1110  FORMAT  (17H  HYPERBOLIC  ORBIT) 

734  NL  «  NL+l 

735  GO  TO  86 

736  85  RP  -  SL/( l.+SQRT ( l.-SL/AX) ) 

737  HP  «  (RP-RE 1/6076.104 

740  HA  •  <AX*$L/RP-RE)/6076.104 

741  WRITE  (6,1106)HP,HA 

742  1106  FORMAT  (21X, 7HPERIGEE,F9.2,3X,6HAP0GEE,F9.2) 

743  NL  •  NLd 

744  IF  (HP.LT.50.)  GO  TO  20 

747  86  IF  (S4.EQ.0.)  GO  TO  70 

752  GO  TO  20 

C 

C  INTEGRATE 

C 

753  60  IF  (S3.EQ.-1.)  GO  TO  70 

756  00  65  1*1,10 

757  CU)  »  YOd)*DT 

760  GO  TO  (61,62,63,64) ,J 

761  61  Yd)  «  Y(  I  )+C(  I  )/2. 

762  0(1)  *  C(I) 

763  IF  (I.EQ.l)  T  «  T+DT/2. 

766  GO  TO  65 

767  62  Y ( I )  *  Y( I ) ♦.29289322* (C (I )-Q( I ) ) 

770  0(1)  *  «58S78644*C( I )♦• 121 32034*0 ( I ) 

771  GO  TO  65 

772  *3  Yd)  •  Y  ( I )  ♦1.7071068*  (C  ( I  )-Q(  I ) ) 

773  Q ( I )  «  3. 4142l36»C(I)-4. 1213203*0(1) 

774  IF  (I.EO.l)  T  ■  T4-0T/2. 

777  GO  TO  65 

10U0  64  Yd)  *  Y(IMC(I)/6.-Q(!)/3. 

1001  65  CONTINUE 

1004  GO  TO  20 

C 

C  DIAGNOSTICS 

C 

1005  70  AVU  *  SCRT(Y(5)«»2+Y(6)*»2^V(7)*»2)-l. 

1006  AVW  *  SQRT(Y(8)»»2+Y(9)«»2*Y(10)*«2)-l. 

1007  SPWU  ■  Y(8)*Y(5)*Y(9)*Y(6)*Y( 10)»V(7) 

1010  WRITE  (6, 11 08) AVU, AVW, SPWU 

1011  1108  FORMAT  ( 1H0, 18HACCURACY  CHECK  -  U,F12.8,3H  W,F12.8,5H  U.W,F12,0) 

1012  CALL  OVERFL  (K) 

1013  GO  TO  (71, 73, 72), K 

'  1014  71  WRITE  (6,1111) 

1015  llll  FORMAT  (9H  OVERFLOW) 

1016  GO  TO  73 

1017  72  WRITE  (6,1112) 

*  . . .  . 
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1020  1112  FORMAT  I 10H  UNDERFLOW I 

1021  73  CALL  OVCHK  IK) 

1022  IF  (K.E9.2)  GO  TO  74 

1025  WRITE  I At  11 13) 

1026  1113  FORMAT  (13H  DIVIDE  CHECK) 

1027  74  LO  «  L 

1030  CALL  CLOCK  (L) 

1031  DL  «  L-LO 

1032  OLS  «  0L/60. 

1033  WRITE  I6t 1114)0LS 

1034  1114  FORMAT  (15H  COMPUTING  TlMEtFS.2»SH  SECONOS) 

1035  GO  TO  1000 
C 

1036  END 


PLUS  Pftrnive,  FPT  (uhu MwTMfc  Panrr  Tts**) 


******  ■ 


•tamatm 


mmxmumiwn  r***?'**u-4#***» 


■T 


APPENDIX  IV 


INPOTS 

A  complete  list  of  i*>p:ts  is  giver  on  the  next  page,  with 
symbols,  definitions  and  uritrf.  They  are  divided  into  two 
groups  corresponding  to  the  two  types  of  data  '.ards,  and  the 
order  is  the  sane  as  that  on  the  data  cards.  The  Indicators  and 
ignition  value  are  defined  as  follows 


Ignition 

Indicator 

Meaning 

Ignition  Value 

0 

Ignition  at  given 
-  Height 

Height  (ft)/1000 

l 

-  Flight  Path  Angle 

Angle  (deg) 

2 

-  Time  after  Previous 

Burnout 

t  iterval  (sec) 

Other 

-  Burnout  of  Previ-ws 

Stage 

Firing 

Indicator 

Meaning 

0 

Direction  of  roll  axis  during  firing 
-  fixed,  tangent  to  absolute  flight  path  angle 

l 

at  ignition 

-  parallel  to  relative  velocity  (aero-stab.) 

2 

-  fixed,  given  elevation  and  azimuth  at  ignition 

3 

-  parallel  to  absolute  velocity 

Other 

-  as  0 

The  layout  of  the  data  cards  is 
the  list  of  inputs.  They  are  of  two 
general  parameters,  and  one  for  each 


shown  on  the  page  following 
types,  one  for  gun  and 
of  the  rocket  stages. 
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All  data  fields  are  five  columns  wide.  Decimal  points  need  not 
be  punched  (and  must  not  be  in  the  first  two  fields).  The 
as aimed  number  of  places  to  the  right  of  the  decimal  is 
indicated.  This  format  will  be  superseded  by  a  decimal  point 
if  punched,  but  the  formal  of  the  input  print-out  will  not  then 
agree.  These  cards  must  be  ordered  according  to  the  sequence 
of  operation  in  the  real  system,  for  each  case;  i.e.,  a  gun 
card  followed  by  a  card  fr  the  first  stage,  a  card  for  the 

second  stage,  etc.  The  number  of  stage  cards  must  be  the  same 

\ 

as  the  number  of  stages  quoted  on  the  gun  card.  The  sets  of 
cards  for  successive  cases  follow  one  another. 

The  integration  interval  is  ideally  one  which  produces  the 
least  errors.  As  the  interval  is  increased,  the  mathematical 
errors  predominate;  as  it  is  decreased,  nund  off  errors  take 
over.  For  most  purposes,  5  seconds  should  be  about  right. 


SYMBOL 

DESCRIPTION 

FORTRAN 

MATH 

M 

Case 

Case  Number 

NS 

ns 

Number  of  Stages 

WT  (NS-tl ) 

wns+l 

Payload  Weight 

(lb) 

SA(NS-hl) 

Sns+1 

Payload  Area 

(ft2) 

DTI 

AtI 

Step  Size  for  Integration 

(sec) 

DTP 

AtP 

Nominal  Output  Interval 

(•ec) 

VR 

VR 

Muzzle  Velocity 

(ft/qec) 

ELG 

Gun  Elevation  from  Horizontal 

(deg) 

an.,  (W',‘  '■*  •• 
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SYMBOL 

DESCRIPTION 

EBggSl 

MATH 

AZG 

A* 

Gun  Asiau  th  froa  North 

(deg) 

ON 

* 

Gun  Longi tude  (Positive  East) 

(deg) 

ATG 

Gun  Latitude  (Poaitiva  North) 

(deg) 

HT 

H. 

Kuala  Heigit  Above  S.L. 

(ft) 

11(H) 

Xln 

Ignition  Indicator 

12  (N) 

X2n 

Thrust  Direction  Indicator 

WT(N) 

Wn 

Total  Weitfit  Before  Ignition 

(lb) 

SA(N) 

Sn 

Cross  8ection  Area 

(ft2) 

EA(H) 

Exhaust  Area 

(ft2) 

DTB(H) 

AtBn 

Burning  Tiaw 

(sec) 

QT(N) 

On 

Ignition  Value 

ELS  (N) 

*en 

Thrust  Elevation  above  Horisontal 

(deg) 

AZS (H) 

A*n 

Thrust  Axiauth  fro*  Track 

(deg) 

FR(N) 

*n 

Mass  Fraction 

SI(H) 

*n 

Specific  Impulse  in  Vacuun 

(sec 

rwmtm 
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INPUT  DATA  FORMAT  AND  DECK  MAKE-UP 


APPENDIX  V 


DECK  HAKE-UP 

The  order  of  the  complete  deck  and  the  details  cf  the 
control  cards  for  use  on  the  IBM  7040  at  McGill  are  shown  on 
the  preceding  page.  The  first  is  for  use  with  a  source  deck 
(as  listed  in  Appendix  III).  Normally,  this  would  be  used 
only  when  making  program  changes;  wL  -is  have  proven 

satisfactory,  the  second  card  should  not  be  punched  NODECK 
as  shown;  an  object  deck  will  then  be  produced.  Hie  second 
shows  the  deck  set-up  for  normal  operation,  with  the  object 
deck.  The  FPT  deck  is  an  object  deck  cf  a  special  routine 
which  allows  unlimited  floating  point  traps;  the  normal 
procedure  is  for  the  monitor  to  terminate  execution  when  the 
recorded  number  of  floating  point  traps  reaches  a  given 
amount. 

On  the  first  card  Mt5S  is  the  estimated  computing  time  in 
minutes  including  compilation,  and  "a"  is  the  estimated  number 
of  pages  of  output.  If  either  of  these  is  exceeded,  the  job 
will  be  terminated.  Ihe  computer  time  may  be  estimated  as 
3.8  +  tf/18  A  tj  (or  tf/10  A  tp,  whichever  is  larger)  seconds 
computing  time  per  case  where  tf  is  the  expected  flight  time 
in  seconds,  plus  compiling  time  of  4  minuses  if  a  source  deck 
is  used.  The  number  of  pages  of  output  may  be  estimated  as 
1  +  tj/55  A  tp  per  case  rounded  to  the  next  higher  integer, 
plus  8  for  the  listing  if  a  source  deck  is  used,  plus  3  for 
bookkeeping. 


APPENDIX  VI 


QUIT VIS 

Ihe  print  out  £  jr  each  case  begins  with  a  complete 
listing  of  all  inputs,  followed  by  the  column  headings 
for  the  periodic  records  which  are  printed  at  intervals 
of  a  tp  and  at  events  Oaut-cb,  ignition,  burnout,  apogee 
and  impac  t ) . 

TIMF  . .  time  in  seconds  from  launch 

HEIGHT  . .  above  sea  level  in  •'autical  miles 

RANG!'  . . .  dt  sea  level  (great  circle)  in  nautical 

mi  ies 

VELOCITY  . .... .  absolute  velocity  in.  tt/sct 

AIRSPEED  . . .  velocity  relative  to  rotating  earth 

in  ft /sec 

PATH  ANGLE  . . .  arct«  »  elevf.tion  of  the  tangent  to 

the  ubsv.iute  flight  path  above  the  local 
\  rizontal  in  degrees 

ELEVATICN  .  angle  of  elevation  of  the  vehicle  roll 

axis  above  the  local  horizontal,  in 
degrees 

AZIMUTH  .  angle  from  the  horizontal  velocity  vector 

to  the  horizontal  projection  of  the 
vehicle  roll  axis,  measured  clockwise  as 
viewed  fnt  above,  in  degrees 

WIND  EL  . .  as  elevation,  for  the  relative  velocity 

in  degrees 

WIND  AZ  . .  as  azimuth,  for  the  relative  velocity, 

in  degrees 

In  addition,  messages  are  produced  at  the  appropriate  times 


to  indicate  igni»l'n,  burnout,  apogee  and  impact.  Except  for 
apigee,  these  are  accompanied  by  latitude  (in  degrees  North) 
and  longitude  (i°  degrees  East). 


2  • 


When  burnout  occur*  at  a  height  greater  than  250,000  ft., 
perigee  and  apogee  heights  are  calculated  and  presented 
In  nautical  miles,  unless  the  orbit  is  hyperbolic,  which 
fact  Is  Indicated. 

Various  error  conditions  are  indicated  by  appropriate 
■Mfe^.4,,  "Attempted  ignition  (of  a)  stage  before  burnout 
(of  the  pi ovioua)  stage"  is  caused  by  incompatible  inputs. 
"Iteration  e",  ’’overflow"  or  ’’divide  check-'  require 

detailed  examicsi ’.on  of  the  program.  "Underflow*  ia  not 
signll leant.  '’Accuracy  check"  shows  the  errors  in  the  unit 
vectors  on  the  U  and  W  axes  and  their  scalar  product  after 
conviction  of  the  case.  These  should  be  of  the  order  ;f  the 

g 

round  off  error,  i.s.  cif  the  order  of  t^/10  A  tj.  An 
arbitrary  limit  of  20G  was  placed  on  the  normal  acceleration, 
and  exceeding  this  figure  results  in  termination  :1  the 
calculation  with  the  appropriate  message, 

A  sample  output  is  shown  on  the  following  pages. 


HARP  TRAJECTORY  PROGRAM  -  R.  H.  MCKEE  CASE  -0  PAGE 
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>In  order  to  assess  the  potential  of  multi  stage  gun-launched  rockets, 
a  study  of  the  vehicle  and  trajectory  parameters  was  undertaken.  A 
digital  computer  program  for  trajectories  was  written  and  was  used  in 
an  experimental  manner  to  approach  optimum  performance  within  various 
sets  of  restricting  assumptions.  The  approach  was  found  to  be 
effective  and  a  useful  orbital  potential  was  demonstrated  with 
reasonable  design  parameters.  ,  \ 
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